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Abstract 

The dynamical likelihood method for analysis of high energy collider events is reformulated. The method 
is to reconstruct the elementary parton state from observed quantities. The basic assumption is that each 
of final state partons occupies a unit phase space. The parton kinematics is statistically reconstructed using 
(a) virtual masses of resonant partons and (b) parton kinematic quantities inferred from observed quantities. 
Generation of (b) is made with the transfer function which is the probability function for parton kinematics 
from a given set of observables. Corresponding to the unit parton phase space, the transfer variable spaces 
are also quantized. The likelihood of the reconstructed state is defined by the Poisson probability for a single 
event with the expected number of event that is the cross section per unit phase space times a luminosity 
factor. Applications of the method to selection of process, parton-observable identification, determinations 
of parton kinematics and dynamical parameters are discussed. 
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1 Introduction 



Dynamical Likelihood Method (DLM) was originally proposed in Ref. [1] and developed in Refs. [2] and [3] as 

a method to determine dynamieal parameters, e.g. masses, decay widths or coupling constants from measured 
quantities (observables). The basic idea is to use the differential cross section ( d.c.s. ) as a theoretical input 
for the event reconstruction. A formulation of the method to use the d.c.s. as a posterior probability was given 
in Ref. [4], and was used in Ref. [5]. 

An alternative formulation of the posterior probability with the d.c.s. is presented in this paper. The 
motivation is as follows. 

In the traditional use of the d.c.s., it is defined per certain kinematic quantities and integrated by other 
kinematic variables. This is to match the detector arrangements devised for individual experiments. The d.c.s. 
there is the prior probability and can be essentially applied to a large number of events. In the general purpose 
collider experiments, however, the detectors are not designed to measure particular quantitites but to make it 
possible to get over-all picture of each event. Hence the event reconstruction on event by event basis is feasible. 
To deal with a single event, we postulate that partons in the final state of hard scattering occupy a unit phase 
space. Thus the d.c.s. is defined per unit phase space of the final partons. 

Leptons, quarks and gauge bosons are called partons in this paper. Parton process described by the Feynman 
diagram, i.e. a hard scatterring process, is called the elementary or parton level process. A process from the 
initial beam state to the observables is called a path. A path that a real event has taken is unknown. The event 
reconstruction is to infer paths which the event could have taken. 

The probability density function ( p.d.f. ) of the first half a path, i.e. from the beam to the parton state, 
is defined by the d.c.s. per unit phase space ( quantum state ) of the final partons in the elementary process. 
This is consistent with the original definition of the d.c.s. 

The p.d.f. for the second half a path, i.e. from the parton state to observables, is the transfer function ( 
T.F. ) that relates the observables with the corresponding variables at the parton level. Corresponding to the 
quantum condition on the final parton state, the transfer variable spaces are also quantized. 

The third quantum condition is about the number of events. The d.c.s. multiplied with a luminosity factor 
can be interpreted as an expectation value of the number of events, and the likelihood ( posterior probability ) 
of the event is defined by the Poisson probability for 1 event with the given expectation value. 

This paper is organized as follows: 

In Sec. 2, the d.c.s. for the final parton state and its use for the likelihood definition arc discussed. The 
parton level likelihood works as the likelihood of the reconstructed parton process as we see in the following 
sections. 

The subject of Sec. 3 is the T.F. A way to obtain the T.F. with fully simulated Monte Carlo events is 
proposed. The detection efficiency associated with measurements and event selection criteria is automatically 
included in the T.F. According to the quantum parton state, the Jacobian scaled transfer variables are intro- 
duced. As a result, the likelihood is essentially Jacobian free. The T.F. in its posterior form is used for parton 
level reconstruction. 

An efficient way of inferring paths is discussed in Sec. 4. Momenta of observed final partons are inferred 
with the posterior T.F. using observed quantities as inputs. Virtual masses squared of resonant primary partons 
are inferred with their propagator factors. 

In Sec. 5, the likelihood for reconstructed multiple parton states is discussed. To regenerate the unknown 
true path of a given event, one makes multiple inference of paths in the event. For the multiple inferences, the 
maximum, the expectation and the multiplicative values of the likelihood are defined. They are used statistically 
to select the process, the parton-observable identification and the kinematic solution of the secondary partons. 
The determination of the dynamical parameters is made by the joint likelihood of all events in the sample. 

Summary of the formulation is given in Sec. 6. 

2 Differential Cross Section and Parton Level Likelihood 

In DLM, the d.c.s. for an inferred parton process is used as a theoretical input to evaluate the likelihood of the 
reconstruction. 
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2.1 Differential cross section for a final peirton state 

We assume that a process is described at the parton level by 

a/A^-b/B^ > ci + • • • + c„ = C, (1) 

where a and h are the initial partons, each representing a quark or an anti-quark or a gluon, in beam particles 
A and B respectively, and c(ci,C2, •.•c„) are the final state partons. States of partons are after the initial- and 
before the final-state radiations. Throughout this paper, particle symbol p also represents its 4-momentum, and 
p its 3-momentum. The final partons are assumed to be on mass-shells, i.e. 3-momenta are enough to define 
their states. Process (1) as a whole, i.e. a set of momenta of all partons, is called parton state in this paper. 

Beam particles A and B are assumed to make a head-on collision along the 2:-axis. Then the hadronic 
cross-section for process (1) is given by 

da = dZadzbdPTfa/A{za, oc)fb/B{zb, a.)fT{PT, a)da{a + b^C; a), (2) 
where da is the parton level cross section, 

daia + b^C;a)= ^^^pl^±l^\M{a + b^C; a)p<i$(/). (3) 
4^y{a ■ 6)2 - mlm^ 

In Eq. (2), symbol a stands for a set of dynamical constants, e.g. masses, decay widths or coupling constants. 
Hereafter, we use symbol a to represent only unknown parameters to be measured. Variables Za = az/\A\ and 
Zb — bz/\B\ arc momentum fractions of a and b in hadrons A and B respectively, and Pt is the total momentum 
of the initial/final system of process (1) in the plane perpendicular to the beam axis. The probability density 
functions ( p-d.f.) for Za, Zb and Pt are denoted by fa/Ai fb/B^ and /t, respectively. Functions fa/A^fb/B and 
/t are effective parton distribution functions for process (1) with the radiation parts removed. In practice, 
these p.d.f.'s are to be obtained by running Monte Carlo event generators. 

In Eq. (3), M. is the matrix element for process (1), and is the differential factor, 

= rr ^^^^ (4) 

1=1 

of the Lorentz invariant n-body phase space element, 

d^n = 5^{a + b-C)d^\i\ (5) 

We call d^n the phase space element ( P.S.E. ) in this paper. 

The d.c.s. for a final state defined by c(ci, • • • , c„) is obtained by integrating Eq. (2) with the initial state 
variables Za, Zb and pr as 

da = I{a,b)\M{a^-b-^C;cx)fd^[{\ (6) 



namely. 



where 



'^'^ I{a,b)\M{a + b-^C;ci)\\ (7) 



= i^La 2 .J a/AiZa,Cx)fb/B{Zb,a)fT{pT,a). (8) 

4\A\\B\y/{a-by -mlml 
The formulation in this paper is to define the likelihood to be proportional to da/d^ii . 

2.2 Resonances ( internal lines ) 

Assume a process where resonance r, which corresponds to an internal line in the Feynman diagram, is produced 
as 

a + b^r + Cj+i + Cj+2 H h Cn, (9) 
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and subsequently decays into channel p as 

r^p■.Cl+C2^ \-Cj. (10) 

Propagator factor If the matrix elements for processes (9) and (10) are given as Mprod and Mdec respectively, 
the matrix element squared for process (1) is factorized as 

\M{a + 6 ^ r + Cj+i + • • • + c„ ^ C; a)\^ = \Mprodfn{sr)\Mdecf, (11) 
where Sr is the virtual mass squared of resonance r, given by 

i 



Sr = 



(E^^)'- (12) 



i=l 



The lowest order approximation for propagator factor Il{sr) is given by the Breit-Wigner form, 

where is the total decay width of r. Improved forms taking the higher order effects into account are discussed 
in Ref [7]. 

In the event reconstruction, if all Cj's are inferred from observables, Sr is given by Eq. (12), while if Sr is 
generated according to the propagator factor, Eq. (13), Eq. (12) gives a constraint to (ci, ■ • • , cj. 

2.3 Parton level likelihood 

Reconstruction of a parton level process by using s^'s according to the propagator factors and Cj's with the 
transfer function will be discussed in the following Sections. In this subsection, wo define the parton level 
likelihood assuming that a path, i.e. a set of parton kinematics P{Pi, ■ ■ ■ , Pn), where N is the total number of 
partons in the process, is given. 

If parton kinematics is given, the differential cross section <; = da/d^n^ can be caluculated. The expected 
number of events with cross section <; is 

IJ, = h<; = li rpr, (14) 

d^)i' 

where h stands for a luminosity factor independent of kinematics of the path. We assume that an event takes 
place according to the Poisson probability of a single event P(l; n), i.e. 

dP{l;iJ,)=P{l;iJ,)dn = P{<;)d<;, (15) 

and define the likelihood of the path by 

4°)=P(l;Ai)=/xexp(-/x), (16) 

where suffix 1 and prefix (0) stand for a single path and the parton level, respectively. The p.d.f. for <; = da/d^n ^ 
is given by 

P(c) = Zii? = h/J- exp{-ii) = ll<; exp{-li<;), (17) 

fromEqs. (14), (15) and (16). 

Events in a data set are mutually independent, hence the number of event distribution in a data set with 
the total number of events Ntot is given as 

f^«„4»' (18) 

Peak value and the normalization of L^^^ Likelihood L^^ takes its maximum value 1/e for = 1, since 

dr^°^ 

^ = (1-M)exp(-M), (19) 
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and it is normalized as 

poo poo 

/ if^dii = IJ, exp{-ij)diJL = 1. (20) 
Ja Jo 

Luminosity factor li The expectation value of ^ = da/d^n^ is 

d(7 \ 1 2 

—jj. )= ? P(<r)rfc =^ fJ,^ exp{-fj,)di^ = -. (21) 

or 

If the integrated luminosity of the data set is Lint, then 

A^tot = iint/^W^. (23) 

assuming the detection efficiency is 1. Then, 



H = 



(24) 



where ctt is the total cross section for the process. In event reconstruction it may be interpreted as a function 
of the mass of particle in search, but since ij„t and Ntot are measured/measurable quantities, we interpret 
them as obscrvablcs which are intrinsic to the given data sample just as the observed kinematic variables are. 
Formulation including the detection efHciency will be discussed in the later section. 



3 Transfer function 

The p.d.f. for the second half a path, i.e. a path from the parton state to the observables, is the transfer 
function ( T.F. ). 



3.1 Observables in collider experiments 

Typical collider detectors have calorimeters and the tracking system. Calorimeters and trackers with tgive 
energies and momenta of particles respectively. A jet is generally identified with a quark (anti-quark) or a 
gluon. An electromagnetic shower associated with or without a track is assigned to an electron or a photon. 
A track passing through calorimeters with a minimum ionizing signal is identified with a muon. We call these 
particles observable partons and their measured quantities observables. 

Electrons, muons and photons These particles are relatively well identified and their momenta are measured 
within the detector resolutions. 

Jets Jets are assigned to quarks or gluons. Measured quantities of jets have uncertainties due to statistical 
nature of parton shower, hadronization, resolution of detectors and jet reconstruction algorithm. 

As for the assignment of jets to partons and the relation between their momenta, we make following com- 
ments: 

(1) There are color flows between these partons, and the fragmentation is not independent among them. 
But the effect can be integrated in the transfer function to be discussed later. 

(2) It is possible that a quark/gluon is observed as two or more jets. If two or more nearby jets are merged 
into one and identified with a single parton successfully, the effect can be remedied. Otherwise, the effect results 
as an inefficiency of the reconstruction and/or a shift in the values of dynamical parameters to be determined. 

(3) In general 4-momentum {Ej,pj) is measured on a jet. For quarks in the final partons, however, we 
assume their pole masses. Hence 3 quantities of a jet are enough to infer the quark 3-momentum. The selection 
of these quantities is not unique but is to be made according to the process, the purpose of analysis and detector 
properties. 
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Missing partons For partons which do not interact with detectors, e.g. neutrinos, the missing transverse 
energy(MET) defined by 

= -4"''^ = -(4"'^ + E /^^) (25) 

is measured, where E^''^^ is the measured total transverse energy flow, E^"'^^ is the sum of the transverse 
energy flow measured by calorimeters, and ^ fxr is the sum of transverse momenta of muons measured by the 
tracking detector. All vectors in Eq. (25) are in the plane perpendicular to the beam-axis. 



3.2 Transfer functions for observable quantities 

For a real event in experiment, the final observables are known, while the parton state in process (1) is unknown. 
The event reconstruction in DLM is to infer such unknown parton state that leads to an observed variable set 
y{yi, • • • , jjNv)- The parton variable set corresponding to y is denoted by x{xi, ■ ■ ■ ,XNy). We call x and y 

transfer variables. 

Prior' timisj'cr function The prior transfer function (T.F.) is a p.d.f. for y when x is given and denoted by 
w{y\x\\ip,cx.), where ip is an integer to specify the process. The probability for {x,y) to be in {dx,dy) is 

dP{x,y\\ip) = w{y\x\\ip, ol) dxdy. (26) 

where 

Nv Nv 

dx= Yl dxm, dy=Y[ dym- (27) 

rn=l m— 1 

If w(y|a;||Q:) > 0, a certain value of y should exist. Hence we require the normalization condition, 

[ w{y\x\\a)dy = 1, (28) 

Jw>0 

for any x with w > 0. 

A typical example of a component of y is the energy of a jet. The T.F. for a jet depends on models of parton- 

shower and fragmentation, the detector response and the jet reconstruction algorithm. Thus it is appropriate 
to derive the T.F. by using Monte Carlo event generators with full simulation, where the momentum of each 
parton and measured quantities associated with it are provided. Events are to be selected with the same criteria 
as applied to real data. 

w(y|a;||Q;) from Monte Carlo events The T.F. is a function of multi-dimensional variables x and y. We 
assume that the T.F. is factorisable as 

Nv 

w{y\x\\ip,cx) = Y[w{ym\x\\ip,a). (29) 

m=l 

To illustrate how to get T.F., we take a simple case where a T.F. for y^ depends on the corresponding 
variable x„i only, and a has a single component a. We denote the T.F. by ■u;(?/|a;||a), abbreviating process 
number ip and variable number m. 

Let Tixy denote the density of generated number of events at {x,y), and that at x. T.F. i«(y|a;||a) is 
defined such that the number of events in {dx, dy) is given by 

dN{x,y) = rixydxdy = rixdx x w{y\x\\a)dy. (30) 

With the integrated luminosity Lint, the number densities are given by 

da , , ,, . f , ^ da 



/da 
rixydy = Lint-^- (31) 



The y dependence of the detection efficiency is included in t/;(y|a;| |a). 
From Eq. (30), the T.F. is given by 



nj{y\x\\a) = ^^yZ:^. (32) 



1=1 



6 



Thus ?7;(7y|.x| ja) is obtained by filling the {x,y) histogram with weight l/rix for each event. Weighting by l/ux 
is to avoid the double counting of the cross section factor which exists in the parton level likelihood. Integrating 
Eq. (32) by y and using Eq. (31), one obtains the normalization condition, Eq. (28). Note that the correction 
for the detection inefficiency associated with measurements and event selection conditions is automatically made 
by deriving the T.F. with the Monte Carlo events. 



3.3 Quantization of the transfer variable space 

( f) 

We consider how the quantum condition A$„ ' = 1 characterizes the transfer variable spaces. This condition 
applies to all final partons, while the transfer variables make sense only for observable partons. Hence we first 
discuss the case of observable final partons. 

Jacobian scaled variables for observable final parton We denote transfer variables of the l-th observable 
final parton by xi and the corresponding observables hy yi { I = 1, ■■■ , N*^^ ). Variable xi is a 3-component 
function of pi, and generally A^*^^ < n, where n is the number of final state partons. For the l-th observable 
final parton, we introduce variables {Xi,Xi) by 



k=l 

where J^i and Jyi are the phase space Jacobian factors, 

= {2tt)^2Ei 



3 ^$(0 

dXi = Y[dXik = -T^dxi = Ja:idxi, (33) 

dYi = T\dYik = -r^dyi = Jyidyi, (34) 
ay I 



diPlx,Plv,Plz) 



d{xii,yi2,xi3) 



(35) 



Jyl — Jxl\x=y (36) 

Obviously, the unit phase volume A$^'^ = 1 corresponds to the unit variable spaces 

AX, = 1, All = l. (37) 

A time-ordered path may be described as follows. A single path specifies a unit phase volume (cell) of 
final parton which one-to-one corresponds to a unit volume (cell) of X;, and picks up that of Yj statistically 

according to the T.F. In other words, by condition A$^'^ = 1, variable spaces Xi and yi are quantized. The 



elements of these spaces become from real (continuous) to countable (discrete) almost-infinite numbers. The 

width of the quantized single path is AX/ = AY; = A$J''' = 1. 

Transfer functions for Jacobian scaled variables We denote T.F. for Jacobian scaled variable [X.Y) by 
VF(l^|X||a). To compare the two T.F.'s, w and W, we again treat a case of a single variable set {x,y) and 
{X,Y). The number of generated event in {dX,dY) is expressed in terms of W^(y|X||a) as 

dN{X,Y)=Lint%W{Y\X\\a)dXdY. (38) 

But dN should be proportional to the outlet path width, i.e. 

dN{X, Y) = JydN{x, y). (39) 
Comparing Eqs. (30), (31), (38) and (39), one gets a scale invariance of T.F., 

W{Y\X\\a) =w{y\x\\a). (40) 
Posterior T.F. Posterior T.F. w(a;|j/||Q!) for a single component set {x,y) is given by 

w{x\y\\a) = , , , (41) 
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and is to be used to infer parton variable x from observable y. The posterior T.F. for Jacobian scaled variables 
W(X|F||a) is obtained by 

W{X\Y\\a) = -. TT./viy|. N .y (42) 

For a given Y, the value of X is to be inferred by the probability, 

dP{X; Y) = W{X\Y) dX. (43) 
Using the scale invariance, Eq. (40), one gets 

dP{X; Y) = ^ w{x\y\\a)dx, (44) 
where Jx is the mean value of Jx defined by 

The domint part of w(a;|?/||a) is symmetric with respect to x — x, where 

X w{x\y\ \a)dx. (46) 



Hence the effect of — Jx is cancelled out in the first order, and J^'^ Jx- In this approximation, 

dP{X- Y) w w{x\y\\a)dx. (47) 

Thus the variable quantization is required only conceptually, and in practice one can use the posterior T.F. 
u;(a;|y||Q;) instead oiW{X\Y\\a). 

Missing final partons The only observable about missing partons are the missing transverse energy. MET. 
The sum of transverse momenta of missing particles, T{Tx,Ty), is inferred with T.F. for MET, 'w{T\ ^t||q:). 
The parton level cross section can be written as 

da da ^ m 

'71(1) = ^TTTy^^^^ " ^*mx)dTx 5{Ty - ^ c*^y)dTy, (48) 

where c*x^iy{i = !,■••, M) are the (a;, y) components of missing partons. ^-functions in Eq.(48) give constraints, 

M M 

'-TOX ~ -^X) C^y = Ty. (49) 

m=l m=l 

if). 

Since the quantization requirement for each missing final parton, A$j^ = 1, is for 3-dimensional variables, the 
requirement is compatible with the 2-dimensional constraint, Eq. (49). To summarize, the MET constraint 
Eq. (49) is free from the quantization condition, and the phase space of each reconstructed parton. whether 
observable or missing, is to be taken as 1. The value of da/d^n^ is evaluated with the transverse momentum 
components determined with constraint (48). 



4 Path Reconstruction 

4.1 Primciry and secondeiry pcirtons in event reconstruction 

DLM is a procedure to reconstruct the parton state, i.e. a set of momenta of all partons, P{Pi, ■ ■ ■ , Pm), 
including resonances and final partons. The parton kinematics is defined in general by giving momenta of n out 
of N partons. We call such n partons the primary partons, denoting them by p(pi , • • • , p„) . Momenta of residual 
partons are determined by the energy-momentum conservation at vertices of the Feynman diagram. We call 
these partons secondary partons. These names are only to specify roles of partons in the event reconstruction. 
The selection of the primary partons is optional, depending on the process and the reconstruction algorithm. 



8 



4.2 Specifications of process, topology and solution 



Given an event with observable set y, there are 3 integers to specify a path. 

Process First, one has to assign physics process ip{ip = 1, • • • , Np) which y came from. 

Topology Some of observed partons in an event cannot be uniquely identified with final partons in the 
elementary process. Examples are the same sign electrons or melons, multiple photons or jets. In the event 
reconstruction, one has to assign some components of observable y to a set of parton species to define variable 
X. We call each set of the parton assignment to y a topology in this paper, and denote the topology number by 
it : it = 1, • • ■ , Nt, where Nt depends on process ip. Variable x and hence the value of the T.F. depend on the 
assumed topology if. 

Solution If a process includes resonance(s), whether daughters are missing or observed, one can infer Sr and 
solve Eq. (12) for momentum component(s) of daughter parton(s). The solutions are sorted by the solution 
number ig : = 1, • • • , AT^, where Ng depends on ip and it- 

4.3 Outline of path reconstruction 

The procedure of a single path reconstruction is summarized below. 

(1) One specifics process ip and infers cx imiformly, 

(2) One specifies topology it, and infers parton kinematics as follows: 

(a) One specifies an appropriate set of n primary partons. If all primary partons are observable, one infers 
their momenta p according to T.F. Totally missing partons are classified to primary partons, and a set of their 
momenta is to be inferred uniformly in their phase space. 

(b) If a resonance is assigned to a primary parton, one infers its invariant mass squared Sr with the propagator 
factor n(,s,.), and determines a secondary parton momentum by Eq. (12). 

Such inferences of variables p and are more efficient than scanning them uniformly. We call such inferences 
importance sampling ( I.S. ). 

4.4 Inference of pctrton momentum from jet 

Quarks and gluons in the final parton state are observed as jets (ji, • • • , jtv ^t)- The parton momentum can be 
inferred from observables of corresponding jet by using T.F. w{yi\xi). In the following, we abbreviate parton/jet 
suffix i. Variable x can be {E,6,^), {ET,r],<p) of the parton, or any other set as long as it is observable and 
determines c uniquely. 

An efficient way of inferring x, a component of x, is to make a variable transformation. 



where m is a normalized uniform random number (n.u.r.n.: < m < 1), and range [xminiXmax] is defined by 
w >Q. Generating u, one can determine x. 

4.5 Inference of missing transverse energy 

We denote the transverse energy flow of the i-th missing parton by ti{tiC0s4>i,tiSin(f)i). The total transverse 
energy of m missing partons is 



A simple example is a case where only one neutrino is involved in the process, where T = vt- T is a parton 
variable to be inferred with the transfer function. 

The choice of the transfer variable set to infer T depends on whether process (1) includes partons going to 
jets or not. Let xt and yr denote such a 2-dim variable set in general. 

(i) Take xt = T and yr =^t, if no jet is involved in the process. 




(50) 




(51) 
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(ii) If jets are involved in the process, the fluctuation of is strongly correlated with that of the jet energy. 
In this case, take 

Njet Njet 

XT=T+^CTj, yT=V^T + ^ETj (52) 

whore j is the jet number, Exj and cxj arc the jet and corresponding parton transverse energy, respectively. 
Parton transverse momenta Ctj 's are independently inferred from jets. 

In both cases, we assume {xtiVt) part of the transfer function can be factored out as 



w{y\x) oc w{yT\xT), 



with a normalization condition. 



VTn 



w{yT\xT)dPyT = 1- 



Inference of Xt is made by a 2-dim n.u.r.n. as 



u 



w{yT\x)dx 



w{yT\x)dx. 



(53) 
(54) 

(55) 



4.6 Inference and use of Sr 

Inference of Sg Wo consider a case where resonance r is selected as a primary parton and Sr is inferred with 



the normalized propagator factor, njv(sr) = NIl{sr), as a p.d.f. for Sr, i.e. 

POO POO 

/ ILN{Sr)dSr = N IL{Sr)dSr = 1, 

Jo Jo 

I r°° MT 

N = 1/ U{Sr)dSr « . 

Multiplying S{sr - ^i)^) dsr{ = 1) to Eq. (7), one gets 



da 



da 

MP 

da 



^{^r - (^Cj)^)rfSr, 



d{u — u)du, 



where 



u{Sr) 

du 



I njv(s)ds = / n(s)ds / / n(s)ds, 

Jo Jo I Jo 

IlN{Sr)dSr, 



(56) 
(57) 

(58) 
(59) 

(60) 
(61) 
(62) 



i=l 



Equation (61) indicates njv(s) is a p.d.f. for Sr- Thus, in the reconstruction, scanning of Sr can be made 
efficiently by generating a n.u.r.n. u {Q < u < 1), and making a variable transformation from u to by Eq. 
(60). 

If there are a total of resonances, \M.{a + 6 — > C)|^ contains propagator factors, and one can choose 

h {h < Ur) resonances as primary partons. 
The d.c.s. in this case is written as 



da 



da 



3p 



i=l 



(63) 
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The values of is inferred with Tl{sr) (r = l,---,h), and h components of daughters, one for each r, are 
determined by solving simultaneous equations, 

Sr-{f24''''f = {r = l,---,h). (64) 

i=l 

( f) 

The vahic; of da/d'^n is to be evaluated using nionic'ntuni components thus determined. 
For multiple resonances, s^-'s can be scanned independently by Eq. (61). 

Sr for observable daughters When daughters of a resonance are all observable, one can evaluate Sr by Eq. 
(12), using Ci's inferred with T.F. and assumed masses of the final partons. 

An alternative way of reconstruction is to infer Sr according to Eq. (61). This is more efficient than scanning 
daughter momenta Cj's independently, because independent scanning of Cj's generally results in off-resonant value 

of Sr- 

An example is process W ^ qq'. We assume that directions of 2 partons are regenerated from those of 2 
jets with their T.F., and ask energies of 2 jets. In this case, one regenerates sw by Eq. (61) and the energy of 
one parton by Eq. (50), then the energy of the other parton is given by solving equation sw = {q + q')^, and 
its T.F. is used as a factor of the likelihood. 

Sr for missing partons We consider a process, where there are m missing partons, ^i, ■ • • , f!m, and 
intermediate partons. The degree of freedom for missing partons is 3m, while mcasiircmcnt of fJx gives two 
constraints. Thus, if rir > 3m — 2, one regenerates s(si, • • • , Sh.){h = 3m — 2) using Eq. (61) and solves Eq. (64) 
for t/. Then all components of t/ are determined. 

If Hr < 3m — 2, the degree of freedom for ^ is 

d = 3m-2-nr > 0, (65) 

and d components of missing partons remain undetermined. 

Examples of d=0 case Examples of 1 and 2 missing particles are given in the following. 

Example 1: Single W ^ lu production associated with/without jets. In this case, m = rir = l(W^), 
hence if we regenerate sw, then h — l,d = 0, and Eqs. (12) and (51) lead to a quadratic equation for Uz- The 
parton kinematics is determined within two-fold ambiguity. 

Example 2: Dilepton channel in tt production, 

tt l+rbbv9. (66) 

For this process, m = 2 {u and V), and rir = 4(t, t, ,W~), hence d = 0, if we regenerate st, Sf, sw+ , s^- by 
propagator factors and T by the transfer function. Six constraints by Eqs. (12) and (51) lead to a bi-quadratic 
equation for and Ep, and the parton kinematics is determined within 4- fold ambiguity [3]. 

Undetermined variables of missing partons(d > 0) There are cases where some components of the parton 
momenta are left undetermined {d > 0): e.g. in search for SUSY particles where many missing particles are 
involved in the process. If a parton momentum contains such component (s), the parton is to be assigned as 
primary, and the component (s) are to be scanned uniformly in the phase space. 

Example 3: Charged Higgs production in tt channel 

ti ^ {bW+){bH-) {bl+v){hT-Vr) (67) 

^ {bl+v){bl-vmDr), (68) 

where I = e or jj. Here, m = 3{i?i , , Vt) , and Ur = 5 {t,t,W^ , H~ ,t~), hence with h = 5, d = 2. To 
determine the kinematics, one regenerates st,St,sw+,SH- ,St with the propagator factors, Tx,Ty with the 
transfer functions, and any 2 (= d) components of neutrino momenta uniformly in the phase space. Then 
neutrino equations are reduced to the case of Example 2. 

5 Likelihood of Reconstructed Paths 

5.1 Likelihood for a single path and multiple paths in an event 

In this subsection, formulas are for each set of (ip, «*,»«), which are abbreviated. 
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5.1.1 Luminosity factor and the likelihood for real events 

To infer a set of single path kinematics, P{Pi, ■ ■ ■ , Pn), for a given event, we use in general virtual mass squared 
of resonances and parton kinematic variables x. For a set of P, we define the likelihood of the path similar to 
lJ"' of Eq. (16). The only modification required for real data is that for event detection efficiency (acceptance). 
Denoting the efficiency by e(^), the luminosity factor h is to be replaced with 

Ntot e(0«^T ^ ' 

and by replacing li with li, the expected number of events has the same form, 

The likelihood for path k in event i is given by 

Li(a,P|y||*,fc)=4°)(/iifc), (71) 

if all components of x are used to define P. If P is defined with unused components of x, x' , the T.F. for these 
components is to be multiplied to likelihood Li{p), namely 

Li(a, P\y\\i, k) = L^°\fiik) w{y'\x'\\cx, i, k) (72) 

5.1.2 Likelihood for multiple paths in an event 

To infer the unknown true path of an event, one makes multiple path reconstructions. Here we discuss three 
kinds of the likelihood for the true path. The advantage of one to the others depends on the process and the 
purpose of analysis. 

Maximum likelihood The M.L.E. of x, P and a in an event, which we denote by x, P and di, are obtained 

by (a) using general purpose minimum search programs for —2ln{Li) or by (2) joint likelihood for multiple 
paths in an event, to be discussed in the following. By the use of P and x, one can define a likelihood for the 
i-th event, as a function of a, Li{a\y, P\\i). 

ExpectaMon value of likelihood The expectation value of the likelihood for cx as obtained by a total of K 
paths for the i-th event is defined by 

Li(a|y||z) = -^Li(a,P|y||z,fc) (73) 
fc=i 

The expectation values of x, Sr, P and a, which we denote by x, Sr, P and ai, are obtained as their means 
weighted by Li{P, Q:|y||i, k). 

Joint likelihood The value of true value of parton kinematics P, Pq, in an event is unknown but common 
to all reconstructed paths in an event, namely, Pq is identified with a parameter set. Reconstructions of P 
can thus be interpreted as pseudo-experiments to determine Pq , where the single path likelihood plays a role of 
p.d.f. for P. Formally, one inserts 5{P — Po)dP (= 1) into the likelihood, interpreting P and Pq as variable 
and parameter sets, respectively. The joint likelihood for K paths. 



Lf\a,P\y\\i) = 



K 



JJii(a,P|y||i,fc) 



.fe=i 



l/K 



(74) 



can be used to get the M.L.E. P by the method of maximum likelihood (m.m.l.)[9]. The likelihood as a function 
of OL with P obtained from the joint likelihood is denoted by L^-^\oL\y, P\\i). 
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5.2 Likelihood for process, topology and solution 



In the preceding subsection, the hkeUhood is for a given set of {ip,it,is) in an event. We consider next the use 
of DLM for selection of these integers. 

Integer likelihood A{ip, it, ig) We denote the Ukelihood for these integers by A{ip, it, ig). The integer hkelihood 
can be normalized as 

Np Nt Ns 

i=i j=i k=i 

Individual likelihoods for ip, it and is are given by 

Ap = ^^A(ip,it,i,), (76) 

it is 

is 

As{ip,it) = A{ip,it,is)- (78) 

Likelihood Aj, is used to discriminate the background against the signal, At (ip) to select topology in a signal- like 
event, and As{ip, it) to choose solution of Eq. (12) for a likely topology in the signal-like event. 

Evaluation of A{ip,ii,is) by DLM The values of A's arc often provided from other information, e.g. b- 
tagging with vertex measurement selects certain processes and topologies. We denote A's from the other 
information by A^°^ 's, and define the integer likelihood as a function of a for the i-th event by 

Ll{ip,it,is,a\y\\i) = A^°\ip,it,is)\i{a.)ip^i^^i^, (79) 

where A, is the likelihood for the multiple inferences in an event as defined in the preceding subsection, 

Ai(a)jp,it,j, = L^^\a.\y,P\\ip,it,is,i), Li{oL\y\\ip,it,is,i) 

or Lf\oL\y,P\\ip,it,is,i)- (80) 

The values of likelihood i*'s defined by Eq. (79) and (80) are functions of a. Thus it is appropriate to take 
their mean value in the search range of a. Denoting their mean values by L^, A's are given by 



A{ip,it,is) = ^ ^ ■ (81) 

[a] 

where sufRx \a] stands for the search region. 

If the search region is wide, the discrimination power for {ip,it,is) is weak. Thus evaluation of A's and 
squeezing the search region of ot are to be alternately iterated. The M.L.E. of a is obtained by using all events 
in data, as we discuss in the next subsection. If values of A's converge after the iterations, statistical selection 
of ip, it and ig can be made. The whole procedure studied with Monte Carlo events can be applied to real data. 



5.3 Mctximum likelihood estimate of ot from multiple events 

The determination of a is to be made by otNev, i-e. M.L.E. from a total of N^y events in the given sample. The 
simplest way is to fit the distribution of a.i for individual events, obtained from Eq. (??), with those of Monte 
Carlo events with known values of a. The minimum of the fit gives djveti [10]. 

Since events are mutually independent, the ocNev search can also be made with the joint likelihood of Nf,y 
events. Namely, otNev is ot that maximizes the joint likelihood, 

Nev 

LNev{ot) = Y[ X] L\{ip,it,is,oc\y\\i), (82) 

»=1 (ip,»ti»s) 

with L\ given by Eqs. (79) and (80). If the selection of {ip,it,ia) is not uniquely made, ajqev determined from 
Eq. (82) is generally shifted from true value olq because of remaining false sets of {ip,it,is) in the sum. This 
deviation is to be corrected by the Monte Carlo simulation. 

As we discussed in the preceding subsection, alternate iterations of the aissgy search and the selection of 
{ip, it, is) are to be made. If the value of otNev converges, it can be used to redetermine P and L\ in each event. 
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6 Summary and Comments 



The dynamical likelihood method ( DLM ) is formulated as a procedure to reconstruct the quantum process. 

General comments on the formulation For a single event reconstruction, we require 3 quantum conditions: 

(1) the d.c.s. is per unit phase space, da/d^n \ (2) the transfer variable spaces are quantized by Jacobian 
scaled variables, (3) the likelihood is defined by the Poisson probability for 1 event. In condition (1) the final 
state density which plays an important role in the traditional use of the d.c.s. is missing. By the Jacobian 
scaled variables in condition (2), the Jacobian factor, i.e. the final state density, is absorved in the quantized 
path, and the use of T.F. with ordinary quantites are justified. The state density is resumed implicitely by 
condition (3), since the number of event distribution, which is the outcome of the traditional form of d.c.s., is 
given by the likelihood of our definition. In short, one can forget the Jacobian factor in the formulation given 
in this paper. Only exception is the totally missing particles, the reconstruction of which should be made per 
unit phase space. The integration by unknown variables is not to be made in this formulation. 

The luminosity factor li is a constant depending on the event detection efficiency. This factor can be obtained 
from the mean value of the d.c.s ( for the reconstructed parton kinematics ) of individual events, the integrated 
luminosity and the total number of candidate events. In this formulation, the absolute value of the likelihood, 
i.e. the coupling constant for the process, and the dynamical parameters are simultaneously determined. 

The formulation is more suitable than the earlier ones, Refs. [5], to analyse events of the collider 

experiments with An detectors. 

Procedure of path reconstruction Given a set of observables of an event, one defines the primary partons 
and infers a path. A path is sorted by the physics process, the parton-observable identification (topology) and 
the solution for the momentum components of the secondary daughter partons. 

Dynamical constants and parton kinematics in a path are inferred by random number generations: (a) 
dynamical parameters uniformly, (b) 3-momenta of observable primary partons according to transfer functions, 
and/or (c) virtual masses of intermediate partons with propagator factors. If there remain undetermined 
momentum components of missing partons, (d) they are to be inferred uniformly in the phase volume of the 
partons. 

Applications Sck'ctions of the process, the topology, and the solution for momentum components of the 
secondary partons, which are specified by integers, are made by the likelihood values for multiple inferences 
in an event. The parton kinematics for each event is given by the M.L.E. or the expectation value in the 
event. Dynamical parameters are given by the M.L.E. from the joint likelihood of all events. Iterations with 
alternate evaluation of the likelihood for the integers and for the continuous variables/parameters are important. 

Finally, we comment on the use of DLM for new particle searches. Most theoretical models of new particles 
provide forms of the d.c.s. that can be used for DLM. In addition, the mass value does not strongly depend 
on details of the parton dynamics, but only on its essential part, i.e. the propagator factor of the particle in 
search. Thus the search for theoretically unpredicted new particles by DLM is also made possible. 

7 Acknowledgment 

The author would like to express his thanks to Drs. Alvin V. ToUestrup, Melvyn E. Shochet and other members 
of the CDF (Collider Detector at Fermilab) collaboration. This paper came out from useful discussions with 
the collaborators. 

Thanks are also to Drs. Y. Hara, I. Ohba, S. Parke and M. Mangano for suggestions and comments during the 
course of development of the DLM formulation. K. Yorita, K. Ebina, R. Tsuchiya, Y. Kusakabe, J. Naganoma, 
S. Toya, K. Ikado and T. Arisawa have made examinations of the method in specific processes, which are crucial 
for the present formulation. 

References 

[1] K. Kondo : J. Phys. Soc. Jpn. 57 (1988) 4126. 
[2] K. Kondo : J. Phys. Soc. Jpn. 60 (1991) 836. 



14 



[3] K. Kondo, T. Chikamatsu, S.H.Kim : J. Phys. Soc. Jpn. 62 (1993) 1177. 
[4] R.H. Dalitz and G.R. Goldstein : Phys.Rev. D45 (1992) 1531. 

[5] Juan Crur Estrada Vigil : liep-ex/0302031. V. M. Abazov et al, Dzero Collaboration: Nature 429 (2004) 
638. 

[6] K. Hagiwara et al. : Phys.Rev. D66, 010001 (2002). 

[7] T. Sjostrand, L. Lonnblad, S. Mrenna, P. Skand: hep-ph/0108264 (2002). 

[8] See, for example, A. Stuart and J. K. Ord: Kendall's Advanced Theory of Statistics, A member of the 
Hodder Headline Group, London. 

[9] Sec, for example, Rcf. [6], or Ref. [8], or H.Cramer: Mathematical Methods of Statistics, Princeton Univ. 

Press, New Jersey (1958). 

[10] F. Abe et al. Phys. Rev. Lett. 73 (1994) 225. F. Abe et al. Phys. Rev. D50 (1994) 2966. 



15 



